Evolution of parton distributions with 
truncated Mellin moments 
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Abstract 

Evolution equations for parton distributions can be approximately diagonalized and 
solved in moment space without assuming any knowledge of the parton distribution 
in the region of small x. The evolution algorithm for truncated moments is simple 
and rapidly converging. Examples of applications are outlined. 



1. Introduction 

AltarcUi-Parisi evolution equations for parton dis- 
tributions can be exactly diagonalized and analyt- 
ically solved by taking a Mellin transform, which 
turns the x-space integro-differential equation into 
a set of decoupled ordinary first order differential 
equations. This procedure, although elegant, suf- 
fers from a practical drawback: Mellin moments are 
integrals over the entire allowed kinematical range 
for the parton momentum fraction, < a; < 1, 
whereas in practice parton distributions are mea- 
sured only in a finite subinterval of that range, say 
xq < a; < xi. Furthermore, solving the equa- 
tion through a Mellin transform somewhat obscures 
the directionality in x-spacc of parton evolution, 
namely the fact that the evolution in Q^ of the par- 
ton distribution at a given value of the momentum 
fraction z depends only on the values of the distri- 
bution for X > z. In this note I present a method 
PI to evolve parton distributions in Mellin space 
without using any assumptions on the small-a; be- 
havior of parton distributions. The method is based 
on the notion of "truncated" moments, defined as 
integrals over a region xq < x < 1. It turns out 
that (non-singlet) truncated moments are coupled 
by evolution via a triangular matrix whose effec- 
tive size can be kept small without loss of accuracy, 
and that can be diagonalized analitycally with lit- 
tle effort; numerical implementation is thus rela- 
tively simple and efficient. The method addresses a 
specific source of uncertainty in the determination 
of parton distributions (the lack of small-x data, 
which are abundant only for F2^p{x, Q^)). This un- 
certainty is apparent in Mellin space, but is present 



also in a;-space evolution codes, which operate by 
choosing a fixed parametrization for the parton dis- 
tribution and subsequently fitting the parameters 
to the data. Any parametrization has a degree of 
rigidity and theoretical prejudice, so that tacit as- 
sumption made about the a; ^ region affect also 
the region where the data are fitted. This kind of 
uncertainty, in principle very difficult to quantify, 
is replaced in the present method by a small, quan- 
tifiable and sistematically reducible uncertainty due 
to the finite size chosen for the triangular matrix of 
couplings between truncated moments. 

2. Truncated moments 

For the sake of simplicity, in the following I 
will consider the nonsinglet quark distribution, 
q{x,Q^). The extension of the method to singlet 
and gluon distributions requires some work, which 
has recently been completed M. Consider the 
evolution equation 
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where the evolution kernel P{z) is a perturbative 
series in as{Q^), and all dependence on Q^ has been 
suppressed. Define the truncated moment of q{x) as 
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Then the evolution equation for truncated moments 
reads 
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where 
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It is easy to show that the A^-th truncated moments 
is coupled by Eq. (^ to all moments with index 
N + K, [K > 0). Furthermore, the strength of this 
coupling decreases rapidly with K. To prove it, one 
expands GN{xo/y) in Taylor series around y = 1. 
This expansion has radius of convergence I—xq, and 
can be truncated, say at order M. Reorganizing 
the polynomial in y — 1 thus obtained to collect 
the different powers of ?/, and performing the y 
integration, one ends up with an equation coupling 
the A^-th moment to the subsequent M moments, 
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with coefficients Cj^ jsi{xq) that are perturbatively 
calculable for arbitrary xq, given the evolution 
kernel P[z). For moderately small xq, say xq < 0.1, 
one gets a very good approximation to the exact 
evolution for reasonably small values of M, say 
M ~ 3 — 5, and furthermore the value of M required 
to keep the desired accuracy decreases with N. This 
can be understood by noting that, for xq = 0, 
GNixo/y) is independent of y, so that for finite 
but small xq, and over most of the y integration 
range, one is Taylor expanding a function which 
is approximately a constant. Furthermore, the 
contribution to the evolution of higher moments 
is suppressed by the fact that they are dominated 
by high values of x, where the parton distribution 
themselves are decreasing as powers of (1 — x). 
Notice that the accuracy of these approximations 
for each truncated moment can always be controlled 
by comparing numerically to the exact evolution, 
given by Eq. (||). 

3. NLO evolution 

Given the approximations discussed above, the 
evolution of a given truncated moment (say the Nq- 
th) can be calculated by solving a linear system of 
differential equations coupling it to moments Nq + 1 
through A^o + M. Schematically, at NLO, 
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a situation which is formally identical to the 
coupled evolution of the gluon and singlet quark 
distributions. Here analytical and numerical 
computations are greatly simplified by the fact that 



the "anomalous dimension" matrix is triangular, 
so that its eigenvalues are simply given by the 
diagonal entries, while the eigenvectors are given by 
a recursion relation, without the need to compute 
any determinants. Denoting by R the matrix 
diagonalizing the LO anomalous dimension C^-^\ 
one can use it to rotate the set of moments qx 
under consideration, as well as the NLO matrix 
C^^^ , according to 
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Then the NLO solution is given by 
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where the apex denotes quantities evaluated at 
the input scale Qo; and 7m = C^mm/^o- 

4. Perspectives 

The method presented in this note has been 
extended to parton distributions with coupled 
evolution 0, and can be straightforwardly 
extended to polarized partons as well. It can 
be applied to derive estimates for all physical 
quantities depending on evolution, unbiased by the 
extrapolation of data to small x. An example 
already studied in |l| is the running coupling itself. 
Other significant example include the first moment 
of the polarized structure function gi, appearing 
in the Bjorken sum rule ||3|], and the behavior of 
the unpolarized gluon distribution at moderate to 
large values of x. Work on these applications is in 
progress. 
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